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Cigarette smoking is a risk factor for type 2 diabetes. Genetic 
variants in the nicotinic acetylcholine receptor (nAChR) genes 
have been associated with smoking phenotypes and are likely to 
influence diabetes. Although each single variant may have only 
a minor effect, the joint contribution of multiple single nucleotide 
polymorphisms (SNPs) to the occurrence of disease may be larger. 
In this study, we conducted a gene-family analysis to investigate 
the joint impact of 61 tag SNPs in 7 nAChRs genes on insulin 
resistance and type 2 diabetes in 3,665 American Indians recruited 
by the Strong Heart Family Study. Results show that although 
multiple SNPs showed marginal individual association with insulin 
resistance and type 2 diabetes, only a few can pass adjustment 
for multiple testing. However, a gene-family analysis considering the 
joint impact of all 61 SNPs reveals significant association of 
the nAChR gene family with both insulin resistance and type 2 diabe- 
tes (bothP < 0.0001), suggesting that genetic variants in the nAChR 
genes jointly contribute to insulin resistance and type 2 diabetes 
among American Indians. The effects of these genetic variants on 
insulin resistance and diabetes are independent of cigarette smoking 
per se. Diabetes 61:1888-1894, 2012 




Type 2 diabetes disproportionately affects American 
Indians. The prevalence of type 2 diabetes is, on 
average, two to four times higher than that in other 
ethnic groups (1). Although lifestyle and environ- 
mental risk factors are believed to be significant contrib- 
utors to the etiology of diabetes, genetic predisposition 
also plays a critical role (2). Recent genome-wide asso- 
ciation studies have identified multiple genetic variants, 
each of which explains only a small proportion of inter- 
individual variability to diabetes risk (3). It is well accepted 
that the etiology of type 2 diabetes involves many genes, 
but a single gene does not cause disease individually; in- 
stead, multiple genes act jointly in the context of complex 
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networks or biological pathways in leading to disease sus- 
ceptibility and disease development (4). A pathway- or 
gene-family approach taking into account the joint effect 
of multiple genetic variants with marginal individual effect 
may capture a large proportion of the associated genetic 
variants and thus should have a higher power than single- 
gene analysis in dissecting the complex genetic etiology of 
diabetes. 

American Indians have the highest prevalence of ciga- 
rette smoking of all U.S. ethnic groups (5). Previous research 
from experimental and human studies has demonstrated that 
cigarette smoking impairs insulin action (6), causes p-cell 
dysfunction, and induces insulin resistance (7), whereas 
smoking cessation improves insulin resistance (8) and re- 
duces diabetes risk (9). Therefore, cigarette smoking is 
believed to be an important risk factor for type 2 diabetes 
(10). However, the association between cigarette smoking 
and type 2 diabetes reported in previous studies was pri- 
marily performed in European or African American pop- 
ulations. It is unclear whether this relationship also holds 
for American Indians. 

Nicotine is the major bioactive component of cigarette 
smoke that leads to insulin resistance (11) and diabetes 
(12). Nicotine acts by binding to nicotinic acetylcholine 
receptors (nAChRs), a superfamily of ligand-gated ion 
channels that are widely present within neuronal and 
nonneuronal cell types (13). Evidence from human and an- 
imal research has documented that genetic polymorphisms 
in nAChRs are associated with nicotine dependence (14) 
and lung cancer (15). However, according to our knowledge, 
no study has yet investigated the potential effect of nAChRs 
genetic variants on diabetes risk. Moreover, existing studies 
focused on single-gene analysis, which is less powerful in 
detecting small genetic effect and cannot capture the joint 
contribution of multiple genes. In this study, we conduct 
a gene-family analysis to examine whether 61 tag single 
nucleotide polymorphisms (SNPs) in 7 nAChRs genes jointly 
contribute to the susceptibility for insulin resistance and 
type 2 diabetes in a large, well-characterized American Indian 
population. 

RESEARCH DESIGN AMD METHODS 

Study population. The Strong Heart Family Study (SHFS) is a multicenter, 
family-based prospective study designed to identify genetic factors for cardio- 
vascular disease, diabetes, and their risk factors in American Indians. A total of 
3,665 tribal members (aged 18 years and older) from 94 multiplex families re- 
siding in Arizona (AZ), North and South Dakota (DK), and Oklahoma (OK) were 
recruited and examined between 2001 and 2003. Detailed descriptions of the 
SHFS protocols for the collection of phenotype data have been described pre- 
viously (16). Briefly, a total of 94 families, including 76 three-generation families 
(26 from AZ, 28 from OK, and 22 from DK) and 18 two-generation families 
(5 from AZ, 8 from OK, and 5 from DK), were included in this analysis. 
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The largest family size is 113 individuals from DK, 61 from OK, and 80 from AZ, 
with an average family size of 38 (37 in AZ, 34 in OK, and 45 in DK). The largest 
sibling size is 9 in DK, 9 in OK, and 10 in AZ. All participants received a per- 
sonal interview and a physical examination. The personal interview used a 
standard questionnaire and was administered by trained study personnel to 
collect data on demographic characteristics, medical history, and lifestyle risk 
factors including smoking, alcohol consumption, diet, and physical activity. 
The physical examination included anthropometric and blood pressure mea- 
surements and an examination of the heart and lungs. Fasting blood samples 
were collected to measure lipids and lipoproteins, insulin, plasma creatinine, 
plasma fibrinogen, and glycosylated hemoglobin, and a 75-g oral glucose tol- 
erance test was performed as described previously (17). Laboratory methods 
were reported previously (16,18). The SHFS protocol was approved by the 
Institutional Review Boards from the Indian Health Service and the partici- 
pating centers. All participants have given informed consent for genetic study 
of cardiovascular disease, diabetes, and associated risk factors. 
Assessments of insulin resistance and diabetes. According to the 1997 
American Diabetes Association criteria (19), diabetes was defined as fasting 
plasma glucose >7.0 mmol/L or post-75-g oral glucose challenge blood glu- 
cose of >11.1 mmol/L or receiving insulin or oral hyperglycemic treatment. 
Impaired fasting glucose was defined as a fasting glucose of 6.1-7.0 mmol/L. 
Fasting glucose <6.1 mmol/L was defined as normal. Insulin resistance was 
assessed using the homeostatic model assessment (HOMA) according to the 
following formula: HOMA = fasting glucose (mg/dL) X insulin (pJJ/mL)/405 (20). 
Measurements of diabetes risk factors. Body weight (kg) and height (cm) 
were measured by trained research staff when participants wore light clothes 
and no shoes. BMI was calculated by dividing weight in kilograms by the square 
of height in meters. Waist circumference was measured at the level midway 
between the lowest rib and the uppermost iliac crest with the subjects standing. 
Hip circumference was measured at the level of widest circumference over 
greater trochanters with the legs close together. Waist/hip ratio (WHR) was 
calculated as waist circumference divided by hip circumference. Cigarette 
smoking was assessed via questionnaire and classified as current smokers, 
former smokers, and never smokers. Current smokers reported smoking s 100 
cigarettes in their lifetime and were currently smoking every day or some 
days. Former smokers are those who had smoked >100 cigarettes but were 
no longer smoking. Never smokers are those who smoked <100 cigarettes 
or never smoked in their lifetime. Pack-years were calculated by multiplying 
the number of packs of cigarettes smoked per day by the number of years 
the person has smoked. For instance, one pack-year was defined as smoking 
one pack (20 manufactured cigarettes) per day for 1 year, or two packs per 
day for half a year, and so on. Alcohol consumption was determined by self- 
reported history of alcohol intake, the type of alcoholic beverages consumed, 
frequency of alcohol consumption, and average quantity consumed per day 
and per week. In this analysis, participants were classified as current drinkers, 
former drinkers, and never drinkers. Each participant received a pedometer, 
instructions for wearing the pedometer, and an activity diary at their clinical 
examination. They were asked to wear the pedometer for 7 consecutive days 
and to record the number of the steps taken daily in the activity diary. Physical 
activity was assessed by the mean number of steps per day calculated by av- 
eraging the total number of steps recorded each day during the 7-day period. 
Tag SNPs selection and genotyping. Sixty-one tag SNPs in 7 nAChRs genes 
(CHRNA3-A6, CHRNB2-B4) from the nAChRs gene family were genotyped in 
3,665 SHFS participants. These genes were consistently reported to be asso- 
ciated with cigarette smoking in previous studies. For tag SNP selection within 
each candidate gene, we used the computer program Haploview 4.2 (21) with 
an r 2 threshold of 0.80 for linkage disequilibrium (LD). The following criteria 
were also considered: minor allele frequency (>5%), SNP location (i.e., coding 
region), and Alumina design scores (quantifying how likely a SNP can be 
genotyped). SNPs that could not be tagged (i.e., singletons) were included as 
long as their design score was >0.15. All genotyping was done at the Texas 
Biomedical Research Institute using the Illumina VeraCode technology (Alumina, 
Inc., San Diego, CA). The average genotyping call rates were 98% for the 61 tag 
SNPs, and sample success rate was >99.5%. 
Statistical analysis 

Single SNP association analysis. We examined the association of each 
individual SNP with insulin resistance or type 2 diabetes, separately, using the 
generalized estimating equation (GEE). The relatedness among family mem- 
bers of the study participants was accounted for using GEE by including family 
identifier as a clustering variable. The model also controlled for other covariates, 
including age, sex, WHR, smoking status (former versus current versus 
never), alcohol intake (current versus former versus never), HDL cholesterol, 
LDL cholesterol, systolic blood pressure, levels of physical activity, plasma 
fibrinogen, and renal function (assessed by estimated glomerular filtration rate 
[eGFR]). Results from GEE were further validated by family-based association 
tests using the computer program FBAT (22). Power of genetic association 
analyses for insulin resistance and type 2 diabetes was estimated using the 



computer program PBAT (23). To avoid population substructure, we first 
stratified the analyses by study center (OK, AZ, and DK) and then combined 
the results from three centers by meta-analysis using the truncated product 
method (TPM; described below). 

Gene-based and gene-family analysis. The association of a candidate gene 
(including all SNPs within the gene) with insulin resistance or type 2 diabetes 
was assessed by combining P values from single SNP association analysis. This 
was done using the TPM, as described below. A gene-family analysis was then 
performed by combing the P value of each gene from gene-based analysis, 
including seven genes in the nAChRs gene family. 

TPM is a P value combination method that is used widely in biostatistics 
(24). Suppose there are N SNPs within a gene, denoted by s 1; s 2 , ■ . ., s N . Let P i 
be the P value obtained from single gene-association analysis for s s in relation 
to insulin resistance or type 2 diabetes. TPM takes the product of the P values 
that do not exceed some prespecified value t, and the test statistic is calcu- 
lated as: 

where I( ) is an indicator function. When the P values are independent, W has 
a known distribution. If the P values are not independent, the distribution of W 
can be estimated using Monte Carlo simulation (24). To calculate W, we need 
to estimate the correlation matrix, for the P values. Define the probits J;, i = 
1, 2, • ■ -N, as fj = $ ~ 1 (p^ where <I> ~ 1 ( • ) represents the inverse of a standard 
normal distribution. Following Demetrescu et al. (25), we assume a constant 
correlation, p, between the probits: cov(tj, tj) = p for i,j = 1, —,N. 
p can be estimated by 

p = max(-^-^, p), 

where 

P = l-^f(^ 2 

1 N 

1 = N^ 

The correlation matrix £ can De estimated using p (26). The choice of trun- 
cation point, t, is somewhat arbitrary, and in our analysis, we used t = 0.1. 
This is based on our previous research demonstrating that TPM is quite robust 
to the assumption of constant correlation. Though different choices of trun- 
cation points have little influence on the results, our simulation studies in- 
dicated that a smaller truncation point may be preferable in terms of type I 
error of the test (26,27). Our previous research has also demonstrated that this 
method has good properties in terms of type I error and power under various 
correlation structures among the P values (26,27). The P value for TPM was 
estimated by 5,000 simulations. To confirm our results by TPM, we also con- 
ducted analyses using other P value combination methods, such as modified 
inverse normal method (28) or Simes' method (29), and obtained similar 
results. In this study, we chose to use TPM because it focuses only on 
P values below some threshold, thus avoiding the offset on the overall effect 
by removing some large P values from the analysis and thus probably gains 
power. 

Other statistical analyses. We first stratified all statistical analyses by 
geographical sites (AZ, OK, and DK), and then conducted meta-analysis to 
combine results from the three study centers using TPM. To examine whether 
cigarette smoking influences the association between genetic variants and 
insulin resistance or type 2 diabetes, we conducted separate analyses to 
compare the results with or without adjustment for smoking status. Multiple 
testing was controlled by Bonferroni correction. That is, for single SNP as- 
sociation analysis, we used the significance level of 0.05/61 = 8.2 X 10~ 4 . The 
significance level for gene-based analysis was set to 0.05/7 = 0.007. In all 
analyses, continuous variables, including insulin resistance, WHR, HDL, LDL, 
level of physical activity, plasma fibrinogen, and eGFR, were logarithmically 
transformed to improve normality. Participants with missing information on 
smoking status (n = 18) or diabetes (n - 39) were excluded from analyses. 
Analyses were done using SAS version 9.2 (SAS Institute Inc., Cary, NC) and 
Matlab 7.10.0.499 (The MathWorks, Inc., Natick, MA). 

RESULTS 

Baseline characteristics of the study participants. 

Table 1 presents the baseline characteristics of the study 
participants by smoking status. Compared with never smok- 
ers, smokers were older (38 vs. 41 years old; P < 0.0001), 
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TABLE 1 

Characteristics of the study participants according to smoking 
status 





Ever smoker* 


Never 


smoker 


P 




(n = 


2,115) 


(n = 


1,532) 


value f 


Age (years) 


41.2 


+ 


15.6 


37.8 


+ 


18.4 


<0.0001 


Male sex (%) 


44.5 




33.9 


i 


<0.0001 


BMI (kg/m 2 ) 
WHR 

Systolic blood pressure 


32.2 


+ 


7.9 


32.3 




7.9 


0.10 
















(mmHg) 


123.3 


+ 


17.0 


121.7 


+ 


17.3 


0.19 


Diastolic blood pressure 


HDL (mg/dL) 


50.4 


+ 


14.4 


51.3 


+ 


14.7 


0.27 


LDL (mg/dL) 


99.4 


+ 


29.7 


96.4 


+ 


28.8 


0.29 


Plasma fibrinogen 
















(mg/dL) 


391.4 


+ 


88.7 


388.1 




92.0 


0.03 


Total cholesterol 
















(mg/dL) 


183.3 


+ 


38.5 


177.1 


+ 


34.7 


0.004 


Fasting glucose (mg/dL) 


115.8 


+ 


52.7 


112.3 


-t- 


53.2 


0.45 


Insulin (|xU/mL) 


18.7 


+ 


20.2 


18.8 


+ 


20.4 


0.49 






1 ° 




2 






Insulin resistance 
















(HOMA) 


5.58 


+ 


7.2 


5.45 


+ 


7.1 


0.65 



Data are mean ± SD unless otherwise indicated. *Former plus cur- 
rent smokers. tP values were obtained by GEE, adjusting for age and 
sex. 



more likely to be males (34 vs. 44%; P < 0.0001), more 
likely to be centrally obese (WHR 0.90 vs. 0.92; P < 0.0001), 
had a higher level of plasma fibrinogen (388 vs. 391 mg/dL), 
and had a higher level of total cholesterol (177 vs. 183 mg/dL; 
P = 0.004). The demographic characteristics of the study 
participants according to study center or diabetes status are 
summarized in Supplementary Tables 1 and 2, respectively 
LD analysis. Information for the 61 tagging SNPs exam- 
ined in this study is shown in Table 2. Pairwise LD patterns 
were examined by Haploview (21) using the definition of 
haplotype blocks proposed by Gabriel et al. (30). The LD 
patterns of the 61 SNPs by chromosome in each study 
center were illustrated in Supplementary Figures 1-3. In 
brief, seven haplotype blocks were identified in partic- 
ipants from OK and DK, and six haplotype blocks were 
found in those from AZ. Two haplotype blocks, one in 
CHRNB2 and one in CHRNB3, were observed in partic- 
ipants from all three centers, but other blocks differ slightly 
according to study center, indicating that genetic architec- 
ture of study participants from the three centers may not be 
exactly the same. We therefore stratified all statistical ge- 
netics analyses by study center. To increase power, we then 
combine results from three centers by meta-analysis. 
Single SNP association analysis. Of the 61 SNPs exam- 
ined, multiple SNPs showed individual association with in- 
sulin resistance or type 2 diabetes. After adjustments for 
age, sex, WHR, smoking, alcohol intake, lipids, blood pres- 
sure, physical activity levels, fibrinogen, and renal function, 
7 SNPs in CHRNA3, 2 SNP in CHRNA4, 7 SNPs in CHRNA5, 
3 SNP in CHRNB2, and 10 SNPs in CHRNB4 were associ- 
ated with insulin resistance, whereas 5 SNPs in CHRNA3, 
2 SNPs in CHRNA4, 7 SNPs in CHRNA5, 3 SNPs in 
CHRNB2, and 9 SNPs in CHRNB4 were associated with type 
2 diabetes. However, only 9 SNPs with insulin resistance and 
10 SNPs with diabetes remained statistically significant after 
correction for multiple testing by Bonferroni. Tables 3 and 4 



TABLE 2 

Information for the 61 tag SNPs in 7 nAChRs candidate genes 
genotyped in this study 



(jrene 


dbSNP identification 

llllllLIJtrl 


\ 1 II (JI l LUaUl I Lc 


A 1 1 cil 


MAF 




1 AC 1 7QA 

rsiuol toi) 


1 c; 
10 


Afhr 


ft (\QA 
U.U04 




rsl loo moU 






n i m 
U.191 




1 OO 1 PiCiQA 


1 c; 
10 


A IC^ 


A 1 OA 

U.19U 




1 on 1 A QOFC 




A IC^ 

A/Cj 


U.U9Z 




rslol /Zoo 


1 c 
10 


A 


n noA 










A 1 AA 
U.1UU 




rso / 4oU / 4 


1 c 
10 


A IC^ 


A AOO 

u.uad 




rso/o/ /o 


1 c 
10 


A IC^ 


A 1 0 1 




rs6495308 


15 


A/G 


0.191 




rsooubo^ 


1 c 

lb 


Aid 


A AOQ 




rs7177514 


15 


C/G 


0.189 


CHRNA4 


rs2236196 


20 


A/G 


0.070 




rs2273504 


20 


A/G 


0.347 




rs3787116 


20 


A/G 


0.066 




rs3787137 


20 


A/G 


0.355 




rs6122429 


20 


A/G 


0.227 














rsl 1637635 


15 


A/G 


0.088 


rsl6969968 15 A/G 0.081 




rsl7483686 


15 


A/T 


0.220 




rsl7486278 




A/C 


0.088 




rs2036527 


15 


A/G 


0.087 




rs514743 




A/T 


0.087 




rs569207 


15 


A/G 


0.190 




rs588765 




A/G 


0.101 




rs615470 


15 


A/G 


0.093 




rs680244 


15 


A/G 


0.102 




rs684513 


15 


C/G 


0.345 




rs8034191 


15 


A/G 


0.090 




rs905739 


15 


A/G 


0.193 




rs951266 


15 


A/G 


0.081 


CHRNA6 


rs2304297 


8 


C/G 


0.316 


CHRNB2 


rs2072658 


1 


A/G 


0.113 




rs2072659 


1 


C/G 


0.024 




rs2072660 


1 


A/G 


0.235 




rs3811450 


1 


A/G 


0.016 


CHRNB3 


rsl0958726 




A/C 


0.325 




rsl3277254 


8 


A/G 


0.325 




rsl3280604 




A/G 


0.325 




rs4950 


8 


A/G 


0.327 




rs4952 




A/G 


0.014 




rs4953 


8 


C/G 


0.014 






rs6474413 


8 


A/G 


0.326 


CHRNB4 


rsl 1633223 


15 


A/G 


0.094 




rsl 1636605 


15 


A/G 


0.255 




rsl2440014 


15 


C/G 


0.254 




rsl2914008 


15 


A/G 


0.006 




rsl316971 


15 


A/G 


0.255 




rsl6970006 


15 


A/G 


0.385 




rsl7487223 


15 


A/G 


0.150 




rsl948 


15 


A/G 


0.090 














rs3813567 


15 


A/G 


0.358 




rs3971872 




A/G 


0.108 




rs7178270 


15 


C/G 


0.094 




rs8023462 




A/G 


0.092 




rs950776 


15 


A/G 


0.083 



MAF, minor allele frequency. 
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show the results of single SNP association with insulin re- 
sistance and type 2 diabetes, respectively, by meta-analyses. 
Results of single SNP association analysis with insulin re- 
sistance and diabetes according to study centers are listed 
in Supplementary Tables 3 and 4, respectively. To confirm 
our results, we also conducted family-based association 
analysis using the computer program FBAT. Compared 
with the results by GEE, the genetic associations with in- 
sulin resistance or diabetes obtained by FBAT are in general 
less significant, and none of these associations survives 
multiple testing correction by Bonferroni (Supplementary 
Tables 5 and 6). 

Gene-based and gene-family analysis. Using P values 
obtained from single SNP association analysis, we conducted 
gene-based and gene-family analysis for both insulin re- 
sistance and type 2 diabetes by meta-analyses to combine 
results from the three study centers. Results for gene-based 
analysis show that, after correction for multiple compar- 
isons, genetic variants in four genes (CHRNA3, CHRNA5, 
CHRNB2, and CHRNB4) jointly contribute to the sus- 
ceptibility for insulin resistance and type 2 diabetes (all 
P values <0.005). Gene-family analysis comprising all 
seven genes demonstrates that the nAChRs gene family 
as a whole is significantly associated with both insulin 
resistance and type 2 diabetes (both P < 0.0001). Results 
for gene-based and gene-family analyses are presented in 
Table 5. Center-specific results for gene-based and gene- 
family analyses are shown in online Supplementary Table 7. 



DISCUSSION 

Using a sample of 3,665 subjects from the SHFS, we con- 
ducted a gene-family analysis to examine the joint impact 
of 61 genetic variants in seven nAChRs candidate genes 
on insulin resistance and diabetes in American Indians. 
We found that, although multiple genetic polymorphisms 
are individually associated with insulin resistance and/or 
type 2 diabetes, only a few SNPs survive adjustments for 
multiple testing. However, a gene-based or a gene-family 
analysis considering the joint contribution of multiple 
SNPs reveals significant associations with both insulin 
resistance and type 2 diabetes. To our best knowledge, this 
is the first study examining the joint contribution of the 
nAChRs gene family to the susceptibility for insulin resis- 
tance and type 2 diabetes in any ethnic group. 

Several aspects of our investigation merit comment. 
First, although a SNP may show no or marginal association 
with diabetes by single gene analysis, the joint impact of 
multiple SNPs within a gene or a gene family on disease 
susceptibility could be large. For example, no SNP in 
CHRNA3 is individually associated with diabetes after 
adjusting for multiple testing, but the joint effect of all 
SNPs within this gene reveals a significant association with 
diabetes (P < 0.0001). This finding is in line with previous 
research demonstrating that a gene-based or gene-family 
approach simultaneously modeling the joint effect of mul- 
tiple SNPs within a gene or a gene family may provide a 
better chance to identify genetic variants that would other- 
wise be missed by single gene association analysis (31). It 
is worth pointing out that the observed genetic associations 
by gene-family analysis may not be driven by the gene 
snowing the most significant association in gene-based 
analysis, because removing the most significant gene from 
gene-family analysis did not change our results. Second, 
genetic variants in the CHRNA5/CHRNA3/CHRNB4 gene 
cluster located on chromosome 15q24 have been consistently 
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TABLE 3 

Association of the 61 SNPs with insulin resistance by meta- 
analysis 



SNP 


Gene 


P value 


1 AC 1 7QA 

rslUol Ml) 


Z" 1 ZTDAT A c> 

LjHKJ\A3 


U.lZUo 


rsl 1637630 


CHRNA3 


0.0166 


1 Oft 1 ftftQ/1 


/ VTD \ T A Q 


U.UooU 


■vr-. 1 Oft 1 A OOC 


/'f TD\T 1 o 

CriKNAo 


U.loOU 


■vr-. 1 O 1 


/ > r r r> \ r 1 Q 

CHKNAo 


ft ftft/IO 


rsio /ooyy 


/ > r r r> \ r 1 o 

L^HKisAo 


ri 1 1 ft/i 
U.11U4 


rsd ( 4oU / 4 


f~iTJ~DAT A O 


ft ftftftO 

{).{){){)& 


™-.O r 7/1 Qft r 7Q 


/ tJTT)\T A G) 

CrlKNAo 


ft nooo 


■».r>FC r 70 r 7 r 7G 

rso/cw /D 


L;HKi\A3 


ft nci o 




S^ZTDAT A O 


ft ftQQ/1 


rsbbUboz 


Z" 1 ZTDAT A O 

LjHKJ\A-3 


A 1 TOQ 

U.l tya 


rs l 1 1 (014 


Z" 1 ZTDAT A O 


A AOQA 

U.UdoU 


rszzdbiyb 


Z" 1 ZTD A 7 A i, 

LjHKJ\A-y 


A A 1 1 A 
U.U114 


*. r -.OO r 7QCA/l 

rszz I dO(J4 


ZTD AT A h 


A ACQO 


rsd its 1 lib 


Z" 1 ZTDAT A h 


A AA^JA 

U.UUoU 


rsd lotidt 


, TTTi AT A /. 

CiiRJ\A4 


A 1 C7rt 

0.15/0 


*.r>£ i oo/i on 

rsblzz4zy 


z" 1 ZTD A 7 A i, 


a 1 ccn 
U.lboU 


rsilbddooo 


z" 1 ZTDAT A z 
LjUKNAo 


A 1 CQO 


rsllbd < bdo 


f^ZTDAT A £T 

LjNKJyAo 


A 1 HOA 

U.lbo4 


rsifaybyybo 


ZTD AT A 

LjtiRNAb 


A 1 A C A 

0.1450 


rsl 7483686 


CHRNA5 


0.0052 


rsl7486278 


CHRNA5 


0.0956 


*- r -.OAQ£CO r 7 

rszudboz / 


f^ZTDAT A £T 


A A 1 OA 
U.Uld4 


w.FC 1 A 1 A O 

rsol4 116 


Z" 1 ZTDAT A £T 




rsobyzu t 


TJTiAl A 


A A 1 QO 
U.UlOZ 


rs5ao 1 bo 


LjtlRNAb 


A 1 f\A O 

0. 104Z 


rsblfc>4 iv 


/TTDM A £T 

LjHKJyAo 


A t\CA£. 

U.Uo4b 


rsbd i id 1 


ZTDAT A 


A A 1 CP. 

U.Ulbb 


rs680244 


CHRNA5 


0.1272 


rssUd4iyi 


f ZTDAT A 

LjHKJ\Ao 


A AOQA 


rsy uo i dy 


ZTDAT A 

LjHRNAb 


A H71 O 

0.0 1 VA 


rsyoizbb 


ZTDAT A 

LjHKJ\ad 


A 1 OOA 

U.ldo4 


rsZd(J4zy i 


Z" 1 ZTDAT A H 

LjtlKNAO 


A OAOC 

U.zUZb 


rszu t Zboo 


/"'ZTDATDGi 

LjHKjyaz 


A 1 1 A Q 

U.l 148 


rszO / zb5y 


LjtlRNtiZ 


<10 


rszu t ZbbU 


f~<ZTDATDG) 


A AOO/t 

u.uyd4 


rs38 11450 


CHRNB2 


0.0348 


rsl0958726 


CHRNB3 


0.2112 


... . 1 QO r 7 r 70C A 

rsldZ l (Zo4 


/"'ZTDATDQ 

LjtiKjytia 


A O 1 A/1 

U.Z1U4 


1 QOQACA/1 

rsidzoUoU4 


/"'ZTDATDQ 


A 1 0 1 O 

U.lolZ 


rs4yt>u 


/"'ZTDATD'D 


A OAAA 
U.ZUUU 


rs4yoZ 


/^ZTDATDO 


A 1 1 TO 
U.ll (Z 


rs4yt>d 


fZTDATDG) 

LjtiKjytis 


A 1 1 CIA 

u.iiy4 


rs4yt>4 


/"'ZTDATDQ 

LjtiKNtia 


A 1 OQO 


™J?J7 A A 1 O 

rso4/44ld 


/"'ZTDATD'D 
LjtlKNtiS 


a 1 nnc 

u.iyyb 


rsll636605 


CHRNB4 


<10 -4 


rsl2440014 


CHRNB4 


<10 -4 


rsl2914008 


CHRNB4 


0.0006 


rsl316971 


CHRNB4 


0.0006 


rsl6970006 


CHRNB4 


0.0086 


rsl7487223 


CHRNB4 


0.0754 


rsl948 


CHRNB4 


0.0832 




rs3813567 


CHRNB4 


0.0002 


rs3971872 


CHRNB4 


0.1102 


rs7178270 


CHRNB4 


0.0268 


rs8023462 


CHRNB4 


0.0782 


rs950776 


CHRNB4 


0.0244 



P values in boldface indicate statistical significance after Bonferroni 
correction. 
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TABLE 4 

Association of the 61 SNPs with type 2 diabetes by meta-analysis 
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P values in boldface indicate statistical significance after Bonferroni 
correction. 



TABLE 5 

Gene-based and gene-family associations for the seven nAChRs 
genes with insulin resistance and type 2 diabetes by meta- 
analysis 
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P values in boldface indicate statistical significance after Bonferroni 
correction. 



associated with nicotine dependence (32,33). Our gene- 
family analysis also reveals a significant association of this 
gene cluster with both insulin resistance and type 2 di- 
abetes, lending further support to previous findings. These 
associations, however, may not be mediated by cigarette 
smoking, because further adjustment for smoking had little 
impact on the observed associations. Moreover, smoking 
status is associated with neither insulin resistance nor di- 
abetes in our study (Table 1), indicating that it is unlikely 
to be a mediator for the association between genetic var- 
iants and insulin resistance or diabetes. Third, previous 
studies repeatedly reported an association of rsl6969968 in 
CHRNA5 with nicotine dependence in European Americans 
or African Americans (15,34). Our analyses, however, did 
not find an association of this SNP with either insulin re- 
sistance or type 2 diabetes, probably due to difference in 
genetic background between American Indians and other 
ethnic groups. It is also possible that this SNP might in- 
fluence diabetes risk through pathways beyond cigarette 
smoking. Fourth, in our study sample, >50% of the study 
participants were smokers. However, the amount of ciga- 
rettes smoked per day is quite low (median = 5). The fre- 
quency of cigarette smoke and the amount of smoke varied 
by study center, with the lowest being AZ, followed by OK 
and DK (Supplementary Table 1). Interestingly, though 
tribal members from AZ had the lowest number of ciga- 
rettes smoked per day, strong genetic associations with 
insulin resistance and type 2 diabetes were observed in 
this study center. This supports the hypothesis that the 
impact of nAChRs genes on the susceptibility to insulin 
resistance or diabetes may not be through cigarette smok- 
ing per se. To examine whether smoking quantity influ- 
ences our study results, we conducted a separate analysis 
by including pack-year instead of smoking status as a 
covariate in the statistical models, but our results remained 
unchanged. 

Smokers tend to be thinner than nonsmokers, and 
smoking cessation results in an increase in body weight 
(35,36). The mechanism underlying the link between cig- 
arette smoking and body weight is unclear, but a recent 
study by Mineur et al. (37) suggests that nicotine may 
stimulate the activity of pro-opiomelanocortin neurons by 
activation of hypothalamic a 3 |34 nAChRs, leading to de- 
creased appetite and body weight. Though smokers tend 
to have lower BMI than never smokers, smokers are more 
likely to have increased central adiposity (38,39), a strong 
determinant for insulin resistance and type 2 diabetes 
(40,41). Furthermore, within a normal BMI, smokers tend 
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to have a greater risk of abdominal fat accumulation com- 
pared with nonsmokers (38,42). In our study sample, 
smokers and never smokers did not exhibit significant 
difference in BMI, but smokers had significantly larger 
WHR than never smokers, indicating that smokers tend to 
be more centrally obese than never smokers. This is con- 
sistent with previous findings (38,42). Our results, how- 
ever, are unlikely to be confounded by obesity because we 
controlled for WHR in all statistical analyses. 

The precise mechanisms through which genetic variants 
in the nAChRs gene family influence diabetes are unclear, 
but because smoking increases inflammation (43) and 
oxidative stress (44), each of which has been strongly im- 
plicated in insulin resistance and diabetes (45,46), it is 
plausible that genetic polymorphisms in nAChRs genes 
may influence diabetes susceptibility through their impact 
on inflammatory and/or oxidative responses to cigarette 
smoking. This hypothesis is supported by our analysis, 
which shows that smokers have a higher level of plasma 
fibrinogen than never smokers. In the current study, how- 
ever, we did not observe a significant difference in the 
prevalence of type 2 diabetes or the level of insulin resis- 
tance between smokers and never smokers, suggesting that 
the impact of nAChR genetic variants on insulin resistance 
and type 2 diabetes may not be due to the direct detrimental 
effect of cigarette smoking per se. Instead, nAChR genetic 
abnormalities may influence diabetes through regulating 
appetite or eating behavior (37,47), thus causing changes 
in body weight or metabolic profiles (48). This hypothesis 
is corroborated by our observation that smokers are more 
centrally obese than never smokers, probably due to over- 
eating. Of course, it is also possible that the nAChRs genetic 
variants could affect diabetes vulnerability through other 
independent yet uncharacterized mechanisms. 

Our study has some limitations. First, though we were 
able to control many of the potential confounders, we can- 
not entirely exclude the possibility of residual confounding 
by other factors such as diet, which is known to be related 
to diabetes (49). However, a previous study investigating 
the relationship between cigarette smoking and glycemia 
indicated that dietary factors may not be a confounder 
(50). Furthermore, it is possible that adjustment for dietary 
factors may be overadjustment, as these could be on the 
causal pathway linking cigarette smoking to diabetes. 
Second, this study used a cross-sectional design, which thus 
precluded causal inference. Finally, our analyses were un- 
dertaken among American Indians, and hence our results 
might not be generalized to other ethnic groups. 

In summary, this study provides initial evidence for a 
joint impact of multiple genetic variants in the nAChRs 
gene family on insulin resistance and type 2 diabetes in 
American Indians. The impact of these genetic variants on 
the susceptibility to diabetes is independent of cigarette 
smoking per se. Our results may provide valuable infor- 
mation for individualized prevention or intervention on 
diabetes in American Indians who suffer from dispropor- 
tionately high prevalence of type 2 diabetes. 
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